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ABSTRACT 

Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. 

The primary results - serving as a follow up to the results of Squire & Bhattacharjee [ arXiv: 1506.04109 (2015)] 

- pertain to the “magnetic shear-current effect” as a viable mechanism to drive large-scale magnetic field 
generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo 
could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion 
disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications 
of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation 
methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo 
are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, 
demonstrating the importance of rotation in determining the relative dominance of each mechanism. 


1. INTRODUCTION 

Understanding the origin and sustenance of astrophysical 
magnetic fields remains an outstanding theoretical challenge. 
Turbulent dynamo, in which chaotic fluid motions act to am¬ 
plify or maintain a magnetic field against dissipation, seems a 
likely explanation, but many questions about such processes 
remain. Interestingly, magnetic fields are generically ob¬ 
served to be correlated over larger scales than the underly¬ 
ing fluid motion, and much of dynamo theory has focused 
on these “large-scale” d ynamos. The well known a-effect 
( [Krause & Radler|['l980| ) may explain such behavior, but re- 


quires some breakage of symmetry in the underlying turbu¬ 
lence (e.g., net fluid helicity). In addition, while the lin¬ 
ear (kinematic) regime of such dynamos may be well un¬ 
derstood, there are s till significant difficulties regarding dy¬ 
namo saturation (see [Brandenburg et al.|2012| and references 


therein) including whether it is even possible for large scale 


quenching ([Kulsrud & Anderson|1992 

Gruzinov & Diamond 

19941 [Bhattacharjee & Yuan|1995[ 

Bo 

Idyrev et al.|2005l Cat- 

taneo & Hughes||2009|). Large sea] 

le velocity shear - ubiqui- 


have some very important role to play. Most obviously, shear 
affects the dynamo through simple stretching of the mean field 
(the Q effect), but a variety of other more subtle ef fects may 
also enhance dynamo action in various ways (e.g., [Vishniac 


& Cho[[2QQH [Blackman & Br andenbui^|2602[ [Brandenburg 

2008} [Tobias & Cattaneol^fTldl ). In addition, shear seems to 


allow the growth of large scale dyn amos without net helicity 
or inhomogeneity in the turbulence ([Brandenburg et al.|2008 


[Yousef et al.[[2QQ8b[ ). Such dynamos may play a fundamen 
tal role in a variety of astrophysical processes where a high 
degree of symmetry is present, for instance, the mid-plane of 
ionized accretion disks. 
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These non-helical shear dynamos have been an object of 
fascination in the dynamo literature for some years. Two 
fundamentally different explanations have been proposed for 
how large-scale fields can be generated without any symmetry 
breaking process. The first - the so-called “shear-current ef¬ 


fect” - is in essence an off-diagonal turbulent resistivity ( [Ra 
gachevskii & Kleeorin[[2003[[2004| . When coupled with the 


shear, even rather small values of this transport coefficient 
can overcome the standard (diagonal) turbulent resistivity and 
cause growth of a mean-field dynamo. The second explana¬ 
tion - the stochastic-a effect - relies on the idea that even if 
the mean ^-coefficients vanish , sufficiently strong fluctuations 
can lead to mean-field growth ([Vishniac & Brandenburg[19'9T 


t’ev|2000|[Heinemann et al.|2011|). This dynamo is not 

L-fieid in the usual sense since it relies on the finite size 


Silant 

mean 

of the system to cause mean-field growth; nonetheless, given 
that the universe is sampling a single realization of turbulence, 
not the ensemble average, such effects could be entirely phys¬ 
ical. At the present time, much of the community appears to 
have converged on the idea that non-helical shear dynamos are 
incoherent in nature; i.e., the stochastic-of effect is more im¬ 
portant than the shear-current effect. The primary reasoning 
is that the crucial transport coefficient required for the shear- 
current effect appears to ha ve the incorrect sign, at least at 
moderate Reynolds numbers ([Radler & Stepanov|2006|[Rudi 


ger & Kitchatinov[[2006[ [Brandenburg et al.[[20b8| ). At the 


same time, given the variety of different, but related , incoher- 
ent dynamo mechanisms that have been considered ([Silant’ev 


2000[ [Heinemann et al.[[2QlH [Mitra & Brandenburg [ [20 H 


Richardson & Proctor|26l2[[Sridhar & Singh|2014| ), it seems 

likely that such effects could be relative ly genericT"^ _ 

Here, following up on a recent paper ([Squire & Bhattachar-[ 
j^[2015b| ), we consider the possibility of large-scale coher¬ 
ent non-helical shear dynamos in the regime of low Reynolds 
numbers. We propose a fundamentally different mechanism 
to those discussed above - that a coherent large-scale mag¬ 
netic field can be excited by small-scale magnetic fluctua¬ 
tions. Why should this be important? In any magnetohydro¬ 
dynamic (MHD) system above low magnetic Reynolds num¬ 
ber, the dynamo at smallest scales in the turbulence grow the 
fastest due to the small-scale dynamo ( [Schekochihin et ak 
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2007| ). Such growth is sufficiently rapid that it always over- 
whelms the large-scale field growth and thus a large-scale 
field must be able to grow on top of both velocity and mag¬ 
netic fiuctuations ([Kulsrud & Anderson|l992} [Boldyrev et al. 


[2005 1 jCattaneo & Hughes||2009| ). This idea is at the heart 
of a-quenching, where the small-scale magnetic fiuctuations 
quench the growth of the large-scale field before it has a 
chance to reach significant amplitude. Our proposal is that 
for non-helical shear-dynamos, the effect of the small-scale 
dynamo is positive, enhancing the large scale dynamo growth 
rate. 

In this paper we focus on understanding such a magnetic 
dynamo in the low Reynolds number regime. In this regime 
the problem becomes substantially simpler, due to the greater 
applicability of quasi-li near approximations and the lack of 
a small-scale dynamo ([Yousef et al.||2008a|). This enables 
the effects of velocity and magnetic fluctuations to be stud¬ 
ied separately (e.g., through driving the induction equation), 
as well as allowing simple calculation of transport coefficients 
and fluctuation statistics. We see that with sufficiently strong 
small-scale magnetic fiuctuations, the character of the ob¬ 
served large-scale dynamo changes, becoming more coher¬ 
ent in time and saturating at higher field strengths. That this 
is a coherent dynamo effect is confirmed through numerical 
evalua tion of the relevant transport co efficients. In a recent 
paper ( [Squire & Bhattacharjee[[MT5b| ), we have considered 
the more relevant case where the magnetic fiuctuations are 
self-consistently excited by the small-scale dynamo at higher 
Reynolds numbers, driving a large-scale dynamo once they 
reach saturation. 

In addition to studying the magnetic dyn amo, we re- 
examine the kinematic dynamos presented in [Yousef et al.j 
( [2008a|b[ ), since it is necessary to understand the intricacies of 
the kinematic dynamo before moving on to the magnetically 
driven c ase. We find that the dynamo seen by [Yousef et ai.[ 
( [2008a[) is indeed a stochast ic-of effect, of the type suggested 
by [Heinemann et ak ( 2Q11[ ), in their non-rotating examples. 
However, anti-cyclonic (e.g., Keplerian) rotation can substan¬ 
tially alter the picture, causing a coherent dynamo to become 
possible by changing the sign of the off-diagonal resistivity. 
This behavior is well explained by the O x /, or R Mler, effect 
( [Krause & Radler|1980t [Moffatt & Proctor} 19 82| ). Although 
not commented on by [Yousef et al.[ ( [2008a[ ), our conclusions 
are entirely compatible with their results, nicely explaining 
the observed trends in growth rates. 

One of our primary motivations in this work has been im¬ 
proving understanding of the dynamo observed in ze ro-net 
flux magnetorotational (MRI) turbulence simulations ([Bran 


denburg et al.[[1995{ [Hawley et al.[[1996t [Lesur & Qgilvie 


2008b| ). Given that such turbulence is simply a shear flow 


with the addition of Keplerian rotation, the results presented 
here should be applicable to some degree. Of course, self- 
sustaining MRI turbulence is highly nonlinear and linear dy¬ 
namo results will be generally inapplicable. Instead, one can 
consider the presence of a large-scale dynamo instability as 
an indication that the turbulence will always be accompanied 
by large scale structures. Given that MRI turbulence is both 
rotating and has strong magnetic fluctuations, it seems rea¬ 
sonable to surmise from the conclusions reached in this pa¬ 
per that a coherent dynamo plays an important role. Fur¬ 
thermore, our recent statistical calculations (see Sec. of 
the nonlinear saturation of unstratified MRI turbulence have 
shown nice agreement with aspects of self-sustained nonlin¬ 
ear simulations, in particular regarding the dependence on Pm 


( [Squire & Bhattacharjee|2015^ . Since a coherent dynamo is 
the only possible mechanism in such calculations, this pro¬ 
vides strong evidence to support the relevance of the magneti¬ 
cally driven shear dynamo to MRI turbulence. Our results re¬ 
garding the MRI dynamo mechanism are qualitatively c onsis¬ 
tent with previous c omputational and analytic studies ( [Lesur 
[& Ogilvie|20Q8a|bj ). For the purposes of understanding MRI 
turbulence, the nonlinear behavior of the dynamo will be im¬ 
portant but we leave this complex topic to future wo rk 
gachevskii & Kleeorin|2004[ [Lesur & Qgilvie|2008a[ ). 


1.1. Outline 

Since results on both the magnetically driven and kine¬ 
matic dynamo are presented, we feel it helpful to provide a 
“roadmap” for paper’s structure. This is intended to outline 
how the central results relate to each other, and convey our 
motivations for structuring the paper as follows. 

As discussed, the most important results of this paper are 
those regarding t he “magnetic shear-current effect ,” which act 
as a follow up to [Squire & Bhattacharjee ( [2015b[ ) in the sim¬ 
pler low Rm regime! However, to be able to convincingly 
interpret results - in particular observations of magnetic dy¬ 
namo in nonlinear simulation - it is necessary to first explore 
the kinematic dynamo, its primary driving mechanisms, and 
its dependence on physical parameters. Thus, we first present 
re sults (Sec. [4]) on the dyna mo mechanism in the simulations 
of [Yousef et al.[ ( [2008a|b[ ), which show that the this kine¬ 
matic dynamo is primarily incoherent (although coherent ef¬ 
fects become important with rotation) and provides a compar¬ 
ison point for later results on the magnetic dynamo. This sec¬ 
tion also acts to illustrate the effectiveness of the quasi-linear 
and statistical simulation methods in disentangling incoherent 
and coherent dynamo mechanisms, and demonstrates that the 
direct measurement of transport coefficients yields results in 
agreement with other methods. 

The magnetic shear-current effect dynamo is then studied in 
Sec.j^ To argue for its existence, we use the same tools as for 
the kinematic dynamo: qualitative examination of the dynamo 
from direct numerical simulation, statistical simulations at the 
same physical parameters as in the kinematic case, and direct 
measurement of transport coefficients. We hope that together 
these methods provide a strong argument for the existence of 
the effect and its potential importance in dynamo theory. 

These sections on the kinematic and magnetic shear dy¬ 
namos are preceded by a theoretical discussion of the dif¬ 
ferent dynamo mechanisms that are possible in this geom¬ 
etry (Sec. 1^, and an explanation of the numerical methods 
(Sec.|3]), including the quasi-linear approximation and statis¬ 
tical simulation methods (CE2). The primary purpose of the 
theoretical discussion is to explain the differences between in¬ 
coherent and coherent dynamos, and what properties might 
be used to distinguish these. A different stochasti c dynamo 
mechanism ( Silant’ev[2000[[^dhar & Singh|2014[ ), based on 
the work of [i^aichnan[ ( |1976| ), is discussed in App. CT where 
we come to the conclusion that this dynamo is unfitly to be 
causing observed field generation due to the effects of off- 
diagonal a fluctuations. We finish the paper in Sec. with 
conclusions, including a detailed comparison with previous 
works, as well as suggestions for future studies. 

Throughout this paper, our nonlinea r simulations will uti - 
lize a similar numerical setup to that of [Yousef et al.[ ( [MQ8a[ ), 
with tall boxes (L^ ^ = Ly) to enhance scale separation, 

and relatively small Reynolds numbers (Re = Rm = 100) to 
avoid the complications of the small-scale dynamo. 
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2. SHEAR DYNAMOS 


In this section we conceptually examine the possibilities 
of incoherent (stochastic alpha), and coherent (shear-current) 
dynamos, arising from nonhelical turbulence in a Cartesian 
shearing box. Specifically, we consider an imposed linear ve¬ 
locity shear, Uq = Sxy, and mean-fields are defined by sim¬ 
ple averaging over the horizontal (x and y)directions. (Note 
that S is defined with a negative sign so as to conform to 
conventions in the literature on astrophysical shear fiows). A 
more comprehensive exploration of po ssible dynamo mecha- 
nisms in this geometry can be found inlMitra & Brandenburg] 

( |20r^ . 


In the conventional way, we start by defining mean and fiuc- 
tuating fields through the relation = Bj-^b = where 
Bt is the full turbulent magnetic field and " is the mean-field 
average (simply a spatial average over v and y). We shall also 
make use of the ensemble mean, denoted (•), which is the av¬ 
erage over an ensemble of realizations at the same physical 
parameters. Averaging the induction equation [see. Eq. 0] 
leads to the standard mean field dynamo equations for the 
mean magnetic field B ( |Moffatt| 197 8} [Krause & Radler|1980| ) 


ture Bi = has the growth rate 


yri 




rjt- 


(3) 


Neglecting by assuming \k^T]xy\ S for all k for 
which scale separation holds, positive dynamo growth 
is possible if -Sr]yx > 0, > k^rj^ The 

maximum growth rate is = \Sriyx\IArit, obtained at 
k = ^j\Sr|yx\|2r|t (if this wavenumber fits in the box). 
Eor a single mode of this dynamo, B^ and By are n out 


of phase, B^ = -k ^J\riyjclS\By, and their phases are con¬ 
stant in time, meaning ReiB^By) = ^{B^BDiByBp. A 

nonzero riyx can arise from the effect of shear, the shear- 
current effect (jRogachevskii & Kleeorin||2003 ), from 
rotation ( Krause & Radler||l9^0t [Radler et al. |2003| ), 
the n X / (6>r Radler) effect, or from a combination of 
both. Since with the shear-current effect, rjyx oc S, the 
maximum growth rate of the coherent dynamo should 
scale as 7 ~ 5^ (this also holds with rotation if kl/S is 
fixed, e.g., Keplerian rotation). 


dtB = Vx(UoxB) + Vx6+ —aB. (1) 


Here S = uxb is the electromotive force, assumed to be of 
the form 5/ = + hijk^j,k + • • • due to scale separation, 

and Rm is the magnetic Reynolds number (inverse normalized 


resistivity). Due to the fact that 
V • 5 = 0 we obtain = 0^ 

components of the 77/,itensorP 

B is 
and 1 

see 

a function only of z, from 
there are onlv 4 non-zero 

Brandenburg & Sokoloff 

(2002|); Radler & Stepanov (p 

JUbf. Expanding Eq. ([JJ) one 


obtains 


dtBx = -ayxd^Bx - (Xyyd^By - rjyxdlBy -r T]yydlBx 

^t^y ~ ^xx^z^x ^xy^z^y ~ hxy^i^x hxx^z^y^ 

( 2 ) 


Stochastic alpha effect: This dynamo arises from the combi¬ 
nation of zero-mean ayy fluctuations and the mean shear 
S. Consider Eq. ^ with rjxy = rjyx = 0, {atj) = 0, 
and again take rjyy = rjxx = For simplicity|^ we 
set axy(t) = ayx(t) = (Xxxif) = 0, assume white noise 
fluctuations in ayy, {ayy {t) ayy {t')) = {a^y)6 (t - t'), and 
again take Bj = Bjoe^^^. One can show usin g standard 
techniques ( [Vishniac & Brandenburg||1997| ) that while 
{Bi (0) decays due to turbulent resistivity, it is possible 
for {BiB*) to grow at the rate 



v / 


where the rfij are defined as the various non-zero components 
of riijk. At this stage, the a or r] are not assumed constant 


m- 


in time, space, or over realizations (i.e., atj ^ 
deed with the mean-field average taken over a finite sized 
domai n they can fluctuate strongly. General symmetry argu- 
ments ( [Radler & Stepano^[2006 [Brandenburg et al.[[2008[ ) 
show that = 0, while there are no such constraints on 


the form of r/ij when effects that break the isotropy of the tur¬ 
bulence are present (e.g., shear, rotation). We shall assume 
that the diagonal components of the resistivity, rjyy and rjxx^ 
are positive, since the scale separation assumptions of mean- 
field theory will presumably become invalid if this is not the 
case. 

The two fundamental dynamo mechanisms we will exam¬ 
ine in this work are: 


Coherent shear dynamo: This dynamo arises primarily 
from the coupling between the off-diagonal resistivity 
T]yx and the shear term -SBx. Specifically, for Eq. ^ 
with atj = 0 and rjyy = rjxx = Tjt (equality of the diago¬ 
nal resistivities is just for simplicity, the dynamo is not 
changed qualitatively by relaxing this), it is straightfor¬ 
ward to show that an eigenmode with the spatial struc- 


^ Specifically rjt^k = 0 since = 0, and rjtjk = 0 if ^ ^ 3 since dxB = 
dyB = 0 . 


Thus, positive dynamo growth is possible if fluctua¬ 
tions in a are sufficiently large. The maximum growth 
rate of this dynamo is y = 0.074 S , obtained 

3X k = ^fS(Dyy|54J]^y^^. Note that in any single re¬ 
alization of this dynamo, as observed in simulation, 
Bx and By will grow approximately exponentially; the 
fact that {Bi) = 0 would only become apparent if a 
large ensemble of simulations were carried out at the 
same physical parameters (with the same initial con¬ 
ditions for the mean field). Importantly, initial condi¬ 
tions must be forgotten over the timescale associated 
with the turbulent resistivity, t ~ (k^T]t)~^ (since {Bi(t)) 
simply decays exponentially), which implies that the 
dynamo cannot have a constant phase as it grows in 
time. Eor a single mode of the dynamo, Bx and By 
are on average 7r/4 out of phase (as for the coherent 

a shear dynamo), Re{BxBy) = 2 “F 2 ^{BxBl){ByBp. 

The stochastic-Qf dynamo will also have a dependence 
on the horizontal domain size, since averaging over 
a larger domain will decrease the size of the fluctua¬ 
tions in a, thus decreasing the magnitude of the growth 

^ When the shear is larger than fluctuations in a, \axy,Q;yx,axx} are each 
subdominant to a^y in their effect on the growth rate, see|Mitra & Branden-| 
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rate. More information, including the effects of other 
nonzero a coefficients and correlations between differ¬ 


ent atj, can be found in Mitra & Brandenburg 
The stochastic alpha dynamo has also been derived 
from the MHD equations directly by quasi-linea rly con¬ 
sidering a collection of forced shear ing waves ( |Heine- 


mann et al.|2Ql 1 [McWilliams |2Q12| ). A fundamentally 


different type of stochastic- a shear dynamo has also 
been proposed a nd studied in |Silant’^ ( |2000| ); [Sridhar 
& Singh (|20141). We explore this further m App. |A| 
where we arrive at the conclusion that the effect is un¬ 
likely to be driving the large scale dynamos studied in 
this manuscript due to the adverse effect of off-diagonal 
a fluctuations. 


Of course, in a real turbulent situation, these two dynamos can 
be mixed, and distinguishing the two may be rattier difficult. 
In particular, the r]ij coefficients discussed for the coherent 
shear dynamo will also fluctuate in time and the mean fields 
will generally be noisy, even if the stochastic alpha effect is 
not the dominant dynamo driver. In this work we will use a 
variety of methods to compare the two in different physical 
situations, from directly calculating transport coefficients, to 
simply observing mean-field temporal evolution. 

It is interesting to note that the growth rate of a stochastic-a 
dynamo can be arbitrarily increased or decreased by chang¬ 
ing the volume of the mean-field average. In particular, an in¬ 
crease in the volume of the average by a factor a must lead to a 
reduction in the magnitude of (a^) by a also, assuming the tur¬ 
bulence in each sub-volume is statistically independent. With 
smaller (a^), a reduction in the dynamo growth rate would 
result. In fact, we see this effect explicitly in the simulations 
presented in Sec.|^by simply doubling the horizontal dimen¬ 
sions of our domain, keeping all other parameters fixed. 

3. EQUATIONS AND NUMERICAL METHOD 

In this section we outline the equations solved, as well as 
outlining our quasi-linear and statistical methods. The fun¬ 
damental equations are the nonlinear magnetohydrodynamic 
equations with a background shear flow Uq = -Sxy and pos¬ 
sible rotation 

dUr dUr 

- Sx^ + (Ut • V) t/r + 2Qz xUt + ^P = 
ot oy 

S Ujxy + -h vW'^Ut (Tu, 

OBj OBj ^ _ 9 

—- SX—— = -S Brxy + V X (t/r X B,) + + (Ti,, 

Ot oy 

V • t/r = 0, V Bt = 0. (5) 


Here Q is a mean rotation of the frame, and v and ff are the nor¬ 
malized viscosity and resistivity respectively. Since all quan¬ 
tities are normalized to one it is convenient to define Re = 1 /v 
and Rm = 1 /17 for the Reynolds and magnetic Reynolds num¬ 
ber. The driving noise {au and ajy) is non-helical and white 
in time, localized in wavenumber around k = Gn with width 
671 / 5 , and is used to generate an homogenous bath of small 
scale velocity and/or magnetic fluctuations Ut and Bj in 
Eq. ^ denote the full turbulent fields (Ut is the velocity not 
including the background shear) - while this notation may 
seem cumbersome, for the remainder of the article we will 
split Bt and Ut into their mean (B = Bt, U = Ut) and fluc¬ 
tuations (m, b). We have deliberately not normalized Eq. 


’ This forcing is of the same form as 


Yousef et al. 


j2008a) . 


with respect to the rotat ion Q as is standard in MRI studies 
( [Balbus & Hawley|1998 ), so as to allow study of shear without 
rotation. Throughout this work we consider initially homoge¬ 
nous turbulence with zero average helicity. We use a Cartesian 
box of dimensions (^Lx, Ly,L^ with periodic boundary condi¬ 
tions in z and y, and shearing periodic boundary conditions in 

V. 

Our primary tool for so lving Eq. 0 is the SNOOPY code 
( |Lesur & Longaretti|2007| ). This solves Eq. ^ with a Eourier 


pseudo-spectral method in the shearing frame, using standard 
methods for dealiasing and remapping. Our standard simula¬ 
tion setup is to seed from random Gaussian initial conditions 
in u and 5 at a very small amplitude and reasonably large 
scales (wavelengths greater than -0.2). The forcing, au (and 
sometimes o"^), causes a small scale turbulent bath of fluctu¬ 
ations, and we study growth of the dynamo on larger sc ales 
than the forcing (i.e., k < 15). As in|Yousef et al.|(|2008a|), the 
separation of scales between mean-fields and fluctuations is 
aided by choosing a box that is very elongated in the z direc¬ 
tion, Lz > Lx,Ly, We study the development of the dynamo 
by numerically averaging B over v and y to obtain the mean 
magnetic fields, B, see Sec.[^ Overall, the numerical setup of 
our non linear runs is nearly identical to that of [Yousef et al. 
( |2008a| ), aside from the addition of forcing in the induction 
equations in some simulations. The Reynolds numbers as de¬ 
fined are with respect to the large scale shear. It is also useful 
to keep in mind more standard definitions of these using the 
small-scale velocity, denoted Rm^ and Re/. Since we use the 
same forcing spectrum throughout this work, these are related 
to Rm and Re through 


Rm/ = 


-Rm = 0.053i/™Rm 




with the similar definition for Re/. Most of the calculations 
presented in this work have Re/ = Rm/ 5. 


3.1. Quasi-linear method and statistical simulation 

Eor certain aspects of this study, we have found it to be 
very useful to study the dynamo using a quasi-linear model 
and statistical simulation in addition to the nonlinear MHD 
equations. Here we outline these methods and the motivation 
behind them. More details can be found in lEarrell & loanno^ 
( |2014| ); [S^quire & Bhattacharj^ ( |2015a| ). 

The basic idea of the quasi-linear model is to split the mean 
field and fluctuations before solving the equations, neglecting 
nonlinearities in the fluctuation equations. The equations are 
thus easily derived by substitution of Ut = U-\-u, Bt = B-\-b 
into Eq. (^, followed by a split of each equation into a mean 
and fluctuating part. This leads to 

(dt - S xdy^ U = -2Qz X U -h S Uxf-i- 

vdlU + (-U ■ V« + b-Vb), 

^d, - S xdy^ B - S B^y + fjd^B + V x{uxb), 

d,B, = d,u, = 0 (6) 
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for the mean-fields, and 

- S xdy^ u = -2Qz XU S u^y + vV^m -Vp- 
(u • Vt/ -rt/ • Vm) -r (5 • -r • WB) -r (7a) 

(dt - S ^^ 3 ;^ b = S b^y + -\-V x{uxB) 

V x(JJ xb) + (Tb, (7b) 

V • M = 0, V • = 0. 

for the fiuctuating fields. Note that only the fluctuations are 
driven by noise. In some cases, it is also convenient to not in¬ 
clude the nonlinear stress feedback on the mean-fields, keep¬ 
ing these static in time by simply setting dtU = dfB = 0 . 
This allows the calculation of the transport coefficients di¬ 
rectly [e.g., Tjyjc, or ayy (t) in Eq. Q] for the chosen form of 
Bx, By, or even U, and is essentially a quasi-linear test field 
method. 


Statistical simulation — Noticing that Eqs. ^ are linear and 
driven by white noise, we can solve for its statistics directly. 
This method h as been termed Stochasti c Structural Stability 
Theory (S3T) ( [Earrell & Ioannou|2Q03|) or the Sec ond Order 
Cumulant Expansion (CE2) ( [Marston et al.||2008| ) - we will 
use the term CE2 in this work. In the context of Eqs. (|^ and 
we form the equation for the second order statistics of u 
and b, 

{uu^) {ub"^) \ 

{bu^) {bb^) r 


C = 


This is 


d,C = m: + + Q, 


( 8 ) 


where {a {t) ait')) = S (t - t')Q (with cr representing both the 
u and b noise) and Ji(U,B) is the linear operator in Eq. 0 ; 
i.e., Eq. 0 is equivalent to 




If we then set / (m, b) = (m, b)^ in Eq. ([^, we can drive the 

mean-fields with the deterministic nonlinear stresses obtained 
through simultaneous solution of Eq. 0 . Note that equating 
f{u,b) with ^/(m,Z>)^ is not valid due to the finite size of 


our system and the fiuctuating nature of horizontal averaged 
quantities (see Sec.[^, we discuss this in more detail below. 

We use a Eourier pseudo-spectral numerical method for 
both direct quasi-linear simulation (DQLS) [Eqs. ( [7]) a nd 
and solving the the CE2 equations [Eqs. ^ and (|^]. The 
codes are written in C++ with MPI paralleliz ation, and use 3/2 
dealiasing with the remapping method of |Lithwi^ p007| ). 
Since one solves for the inhomogenous fiuctuation statistics 
in z, the CE2 code requires a grid of size x Ny x (mN^f 
(where m is the number of variables), and CE2 calculations 
can be relatively expensive. The c odes have been verified 
and test ed in a variety of ways; see [Squire & Bhattacharjee| 
( |2015aj ) for more information. 


Discussion of the quasi-linear method — The quasi-linear 
model involves a rather drastic approximation to the full non¬ 
linear equations. What do we gain by studying such a system? 
Generally, such models have allowed a simpler interpretation 
and study of large scale structure growth in turbulence, and 
been rather useful in a variety of geophysi cal, plasma, and 
fiuid dynamics problems (see, for example, [Earrell & loan- 


nou|20Q9l|2012t|Tbbias & Marston|2Q13l[Parker & Kroinme^ 


2013[ ). In previous work ( Squire & Bhattacharjee|20I5a[ ) we 

have found surprisingly good agreement between saturated 
states of CE2 for the MRI system and nonlinear simulation, 
in particularly a strong scaling with magnetic Prandtl number. 

In the context of the work presented here, the methods pro¬ 
vide a simple way to calculate transport coefficients by fix¬ 
ing the fields in both the magnetically driven and kinematic 
cases, followed by an unambiguous check that a mean-field 
dynamo can be observed at the same physical parameters. 
While versions of the (nonlinear) test-field meth od exist that 
explicitly take into acc ount magnetic fiuctuations ( [Rheinhar^ 


& Brandenburg [[MT0| ), these are relatively complicated and 


in the early stages of development. However, the most im¬ 
portant benefit of the quasi-linear methods is afforded by the 
comparison between CE2 and DQLS. This provides an un¬ 
ambiguous test of whether the dynamo is coherent or inco¬ 
herent, since statistical averages are inserted directly into the 
CE2 mean-field equations and an incoherent dynamo is not 
possible. Thus, if similar results are observed between CE2 
and DQLS, we can be sure that the dynamo arises through 
Pij transport coefficients. Another interesting aspect of CE2 
is that long periods of exponential growth in the mean-field 
can be observed, even when strong small-scale magnetic fiuc¬ 
tuations are present (e.g., due to magnetic driving). This is 
in contrast to DNS or DQLS, where it is generically diffi¬ 
cult to observe exponential dynamo growth in the presence 
of strong magnetic fiuctuations, since the finite size of the do¬ 
main causes the mean-field to come into near equipartition 
with the fiuctuations almost instantaneously. 

Einally, we note that CE2 calculations in fixed mean fields 
are in essen ce the same calculation a s the semi-analytic result s 
presented in [Sridhar & Singh[ ( |2010[ ); [Smgh & Sridhar[ ( [201l} . 
Our results agree with their findings in the 
with only momentum equation forcing. 


; rotation-less case 


4. KINEMATIC DYNAMO 

Before exploring the dynamo with magnetic fiuctuations it 
is import ant to fully understand the kinematic dynamos pre¬ 
sented in [Yousef et al.[ ( [2()08a[ ). With this aim, we have re¬ 
produced many of their simulations across a variety of 5, Q 
and L^, to better understand the fundamental dynamo mecha¬ 
nisms. We present the most relevant of these results here. 

Eor the kinematic dynamo we drive only the momentum 
equation in Eq. ^ (i.e., (Tt = 0), and at these Reynolds num¬ 
bers the small scale field arises purely from tangling of the 
mean-field by velocity, V x (m x B), an effect which is quite 
distin ct from the small-scale dynamo ( [Schekochihin et aT 
2007| ). In both the rotating and non-rotating cases, we see a 
mean field dynamo above some threshold in and (cr^) i.e., 
the dynamo is only excited in a sufficiently tall box if driven 
hard enough. Given the scaling of the growth rates in Eqs. 
andj^and the fact that rjt is the sum of a turbulent and physical 
resistivity, this behavior is expected for both incoherent and 
coherent dynamos]^ Our main finding is that the non-rotating 
case is a stoc hastic alpha dynamo (es sentially that explored 
analytically in [Heinemann et al.|2011[ but including non-zero 
kf) but that rotation qualitatively changes the mechanism, de¬ 
creasing T]yx to negative values and causing the dynamo to be- 

^ One might expect the dynamo to disappear again if <a-^> is increased 
further, due to the increase in rjt causing the dynamo to become stable. This 
behavior is seen in the quasi-linear case, but it seems that at these parameters 
in the nonlinear runs, a small-scale dynamo is excited before this occurs. 
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Figure 1. (a-c) Illustration of By (z, t) from non-rotating turbulence with 
S = 2, Lz = 16, Urms = 0-8, for DNS, DQLS, and CE2 from (a)-(c). In 
the two direct runs (a) and (b), we remove the exponential growth (i.e., plot 
e~^^By (z, t) where y is the measured growth rate) so that the full time evo¬ 
lution can be observed, (d) Growth in time of the mean field for the nonlin¬ 
ear equations (solid,blue), quasi-linear DNS (dashed,red) and CE2 (dash-dot, 
yellow), each at the same physical parameters as in (a). While both nonlinear 
and quasi-linear DNS exhibit a positive mean-field dynamo, the CE2 calcu¬ 
lation does not, illustrating the dynamo must be incoherent. The dotted black 
line shows the energy of u fluctuations. 


come more coherent. 

In both the non-rotating and rotating cases we run DNS, 
DQLS, and CE2 calculations at identical parameters. The pur¬ 
pose of this comparison is primarily to illustrate the difference 
between CE2 and DQLS (due to the incoherent mean-field dy¬ 
namo), while showing that the nonlinear case exhibits a qual¬ 
itatively similar dynamo to DQLS. Although the spatiotem- 
poral evolution of the mean-field is similar in each case, we 
shall see that the DNS and DQLS runs exhibit slightly differ¬ 
ent growth rates. This can be attributed to inaccuracies in the 
quasi-linear approximation at these Reynolds numbers. 


NON-ROTATING DYNAMO Eig. [^illustrates the growth of the 
nonrotating dynamo using DNS, DQ LS and CE2, at S = 2 
and Lz = 16, = Ly = 1. As in [Yousef et al.| (|2008a|), 

we use a resolution (32, 32,512) for DJNS and DQLS, but use 
(32,32,256) for the CE2 run since the algorithm scales with 
so is quite computationally expensive H 
Eirstly, it is worth noting that the mean-field, as plotted in 


Figure 2. Same as Fig.^but for Keplerian rotating turbulence, Urms = 0-75 
(velocity flucutations are suppressed slightly by the rotation). In contrast to 
Fig. IT] the CE2 calculation also shows a growing dynamo, albeit at a much 
smaller growth rate, illustrating that the dynamo is partially coherent. 


Eig.[2a-b), is truly a “large-scale” dynamo. We can estimate 
the wavenumber of By as approximately 3 x 271/L^ ~ 1.2, far 
smaller than the forcing scale, kf = Stt. Next, let us com¬ 
pare the CE2 with the the nonlinear and quasi-linear DNS. It 
is evident that the dynamo in this case is purely incoherent - 
while slow mean-field growth is observed in DNS and DQLS, 
the magnetic field simply decays in the CE2 simulation in ex¬ 
actly the way that would be expected due to a positive rjyx 
coefficient. It is also worth noting the qualitative appearance 
of the mean-fields, which appear to wander randomly, as ex¬ 
pected due to a stochastic-a effect. A final piece of evidence 
for the incoherency of this non-rotating dynamo comes from 
doubling the box size in the v and y dimensions, keeping all 
other parameters fixec|j(not shown). This causes the growth 
rate of the mean-field dynamo to change from y = 0.062 [for 
the dynamo in [Ja)] to being almost stable, y = 0.0096, and 
since a coherent dynamo should be mostly unaffected by such 
a change (unless the added wavenumbers significantly affect 
the transport coefficients), this constitutes a simple check of 
the dynamo’s incoherency without using of the quasi-linear 
approximation. 


^ We have verified that identical results are obtained at half this resolution 
and are confident that = 256 is sufficient to resolve all important scales. 


^ We would like to thank A. Schekochichin for suggesting this numerical 
experiment. 
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Figures. Spatiotemporal evolution of log^o(l^);(z,0|) at = 1, L = 32, 
for (a) nonrotating case {urms ~ 0-47), (b) Keplerian rotation (urms ~ 0-43, 
again velocity fluctuations are slightly suppressed by rotation). The log color 
scale is chosen so as to easily see the mean-fleld phase evolution. The differ¬ 
ence in the two dynamos is evident from the evolution of the phase of By as 
the dynamo grows. While in the non-rotating case the phase wanders some¬ 
what randomly, as is characteristic of an incoherent dynamo mechanism (see 
Sec. I^, we see a relatively constant phase of By in the case with rotation. 
Note mso the faster growth rate of the rotating dynamo. 
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Figure 4. Growth rate of Brms as a function of 2Q./S , for fixed shear S = 2, 
and velocity driving ^ 1 (except for the point Q = 0, for which 
u^ms ~ The shaded region shows where the flow is hydrodynami- 

cally unstable (neglecting dissipation), and the dashed vertical line shows 
the SOCA prediction for where the coherent dynamo growth rate vanishes. 
Of course, due to the strong stochastic alpha effect, the dynamo can still grow 
even when the predicted coherent growth rate is zero or negative. The dotted 
line is an approximate fit of predicted growth rate, Eq. to the data, using 
J]t = i^rmsPkf = 0.018, rjyx = 0.0007 x (20 - S), and = 6.2 x 10“^. 
Error bars are estimated by fitting the growth rate to half of the time-series 
data for each run. 


ROTATING DYNAMO In Fig.|^ we illustrate the same calcula¬ 
tions as Fig.[2 but with a Keplerian Coriolis force [Q = 2/35 
in Eq. ©] a®ed. While the dynamo in the quasi-linear and 
nonlinear direct simulations are similar (with a slightly higher 
growth rate) to the non-rotating case, the CE2 dynamo is 
markedly different, exhibiting mean-field growth. This illus¬ 
trates that adding net rotation to the system enabled a coherent 
dynamo, which can be understood as arising from a change in 
sign of rfyx (see also Figs. [5|[6| below). Th is effect is simply 
the well-known Radler, or x J, effect ( [Krause & Radler] 
|1980t|Moffatt & Proctor|1982|). This idea seems to have been 
missed in [Yousef et al.| ( |2008a| , who state “There does not 
appear to be much difference, qualitative or quantitative, be¬ 
tween the rotating an d nonrotating cases.” The finding agrees 
with analytic results ( [Squire & Bhattacharjee[[2Q15c[ ), which 
show that the contributions to rfy^ from rotation and the shear 
have identical forms, and together give riyx ^ S - 2Q. For 
Keplerian rotation this is slightly negative, leading to the pos¬ 
sibility of coherent dynamo growth. Finally, we have again 
doubled the horizontal dimensions of the box for this rotat¬ 
ing case (not shown), which causes the dynamo growth rate 
to drop from y = 0.067 [in Fig.j^a)] to y = 0.041. A com¬ 
parison with the results in the previous paragraph (y = 0.062 
and y = 0.0096 in the narrow and wide boxes respectively) 
shows that although in the narrow box (L^ = Ly = 1) rotation 
causes only a minor difference to the growth rate (because the 
stochastic-a effect significantly overwhelms the coherent dy¬ 
namo), in the wider box (where fluctuations in a have been 
significantly reduced), the difference in growth rates is much 
more substantial. This behavior is consistent with the rotating 
dynamo being driven through a combination of stochastic-a 
and coherent effects [see Eq.j^below], with the coherent effect 
being mostly unmodified by me change in box dimensions. 

While the By(z, t) evolution pictured in Figs.j^a) and[^a) 
looks qualitatively rather similar between the rotating and 
non-rotating runs, this is not always the case. In Fig. we 
compare spatiotemporal evolutions of By{z, 0. in a longermox 
(L^ = 32) with less driving noise, which causes a lower 
growth rate and a decrease in the relative important of the 
stochastic-Q' effect compared to the coherent dynamo. As 
is evident, the two dynamos are qualitatively different, with 


the phase of By wandering quasi-randomly in the non-rotating 
case, while in the rotating case it is approximately constant 
in time. This constant phase is not consistent with a dynamo 
driven purely by the stochastic a effect (see Sec.®. 

From Figs. and we thus interpret the Ke^rian rotat¬ 
ing shear dynamo around these parameters as being driven by 
both incoherent and coherent mechanisms. This interpreta- 
tion is entirely consiste nt with all numerical results given in 
Yousef et al.[ p008a|b[ ). In particular, their Fig. 5 illustrates 
that the addition of rotation enhances the growth of the dy¬ 
namo in all cases. Furthermore, while y ~ 5 for the non¬ 
rotating dynamo, with rotation it is evident that the growth 
of y is somewhat faster than linear in S . Since one expects 
y ~ 5 ^ for a coherent dynamo (since rjyx itself must scale lin¬ 
early with S for small 5), their observed trends are consistent 
with the dynamo being driven through a mix of incoherent and 
coherent mechanisms. We note that a consideration of wider 
boxes would increase the importance of the coherent effect in 
comparison to the incoherent effect, widening the difference 
between rotating and non-rotating dynamos. 

VARYING THE ROTATION As one final test of the importance 
of net rotation in this system we have run a series of simula¬ 
tions, increasing the rotation from Q = -1 (cyclonic rotation) 
to Q = 4 (anticyclonic rotation). Results from this series of 
simulations are illustrated in Fig. As expected, we see a 
substantial increase in dynamo growth rate as the rotation be¬ 
comes anticyclonic, in broad agreement with the SOCA pre¬ 
diction Tjyx oc S - 2EI. Due to the presence of the stochastic 
a effect, one would not expect a linear scaling of y with Q. 
Instead, the growth rate (including an rjy^ and fluctuating ayy) 
is the most positive root o^ 

-4{al) k^S 2 - 4r,y,k^S^ + f = 0, (9) 

where ^ = 27]tk^ + y ( |Mitra & Brandenburg||2012| . We plot 
a fit of Eq. (^ to the data in Fig. H [with rjt = Urmsl'^kf, 
rjyx = %jco(20 - 5)1, illustrating gocW agreement away from 

^ In reality, one will also see a change in growth rate due to rjxy and fluctu¬ 
ations in the other a coefficients, but these effects seem minor and are ignored 
here. 
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Figure 5. Transport coefficients for the kinematic non-rotating dynamo rixx 
(solid line and circle markers, blue) and rjyx (dashed line and square markers, 
black) as a function of Re (at Pm=l), for =2. The curves show the quasi- 
linear results, calculated using CE2, while the markers show the nonlinear 
test-field calculations with error bars (see text). As is common, coefficients 
are normalized by the “high-conductivity” SOCA turbulent resistivity t/q = 
UrmsIQkf). Across all simulations, the absolute level of the forcing (i.e., cr„) 
is kept constant at the same level as Fig. [T] which means that the lower Re 
simulations have somewhat lower Urms- 

the instability boundaries (Q = 0 and 0 = 5/2). Close to the 
boundary, it seems that some other nonlinear effect may be 
important, increasing the growth rate on the cyclonic side and 
decreasing it on the anticyclonic side. 


4.1. Direct calculation of transport coefficients 

To validate and quantify the conclusions discussed above, 
in this section we directly calculate the transport coefficients, 
comparing results from CE2p^ an d the test-field method 
([Brandenburg & Subramanian|20051) (implemented within the 
framework of the SNOOPY code). The CE2 calculations are 
carried out by fixing the mean fields. By = Byo cos(27rz/L^), 
and driving linear fluctuations to calculate their statistics and 
thus the transport coefficients. These calculations, since they 
are quasi-linear in the shear ing frame, are fundamentally th e 
same as those presented in [Sridhar & Subramanian| ( |2009| ); 


Sridhar & Sir^ ^ ( |201Q| ) ; |Singh & Sridhar] ( |201 1| ) for the non 
rotating casej^ 

Test-fiel d method calculations are carried out in the stan¬ 
dard way ( [Brandenburg et al.]|2008| ) by solving the momen¬ 
tum equation with no Lorentz force, using this velocity field 
to drive a small-scale magnetic induction equation 

dtb‘1 = Vx[m xB‘ 1+ Uxb‘i+ iuxb‘i - u x b‘i)\+fiAb (10) 


for a set of test-fields and specified mean-fields (we 
chose a sinusoidal form for B^). There is no small-scale dy¬ 
namo at these parameters, which simplifies the calculation 
since arises purely due to the presence of B^. Calculations 
are run from ^ = 0 ^ 1000, with the error in the transport 
coefficients estimated by dividing the time-series into N seg - 
ments and calculating the standard deviation of the meanp^ 


DQLS gives identical results to CE2, albeit with errors due to the ran¬ 
dom noise. 

The only substantial difference is the forcin g - a singular forcing 


(Tu ~ d (A: - A:/) is used in Singh & Sridhar 
forcing as detailed in Sec.p] 
riviaTt 


(2011), while we use the same 


Ut course, our calculation is numerical rather 
than analytic and it is triviaTto add the effects of rotation (this is difficult an¬ 
alytically a lthough perturbative methods may be feasible, see |Leprovost &| 
|Kim|2008). 

We generally take N = 100, but results are quite insensitive to this 



Results are illustrated in Figs. and We see that in 
both cases the quasi-linear and nonlinear coefficients agree 
at lower Rm as expected, diverging somewhat past Rm >70. 
In agreement with our conclusions from simulations earlier in 
the section, py^ > 0 in the non-rotating case, while py^ < 0 
with rotation, showing that a coherent dynamo is possible at 
sufficiently small It is also worth noting that the magni¬ 
tude of Pyx is less in the rotating case, as known from SOCA 
calculations ( jRadler & Stepanov|2QQ6[ ^quire & Bhattachar " 


je^|2Q15c'j ). For Rm = 100, as used in Figs. to there 
are some differences between quasi-linear and nonlinear re¬ 
sults due to inaccuracies in the quasi-linear approximation, 
which explains the discrepancy in dynamo growth rates ob- 
servecj^in Fig.[^ Interestingly, given the controversies sur- 
roundmg the kine matic “shear-current” effect ( jRogachevskiij 
& Kleeorin|2003|), nonlinear corrections appear to be particu¬ 


larly important for pyx without rotation (this coefficient shows 
the largest discrepancy between the nonlinear and quasi-linear 
calculations). 

In addition to the results for pyx and pxx shown, we have also 
calculated pxy and pyy by setting Bx = Bxq cos (kiz). By = 0. 
We find that pxx = Pyy to a high degree of accuracy, while 
Pxy, which is positive, is mostly unaffected by rotation. Its 
magnitude (compared to the other p) depends strongly on the 
shear and Reynolds number. Due to the dominance of the 
shear, such an pxy has little effect on the growth rate, even 
though its magnitude is larger than that of pyx. In addition 
to the results illustrated and discussed above, we have also 
verified the expected linear dependence of pyx on S at low 
Rm and confirmed that the transport coefficients change very 
little with Lz over the range = 1 ^ 8. 


5. MAGNETICAFFY DRIVEN DYNAMO 

Having now broadl y understood the shear dynamos of 
[Yousef et al.| ( |2008b|aj ), we examine the effect of small scale 
magnetic fluctuations. Before presenting numerical results, 
it is helpful to explain in more detail exactly what we is 
meant by a magnetically driven li near dynamo. Similar ideas 
have been consid ere d before (e.g.,[Rldler et al.|2003t[Park & 


Blackman|2012[ ), see |Rheinhardt & Brandenburg[ ( [2010[ ) toFi 


particularly thorough analysis. 

As is obvious from Eq. ([T]), an unstable dynamo requires in¬ 


choice. 

It seems that for Fig.fq the lower values of rjxx and rjyx in comparison 
to the quasi-linear runs can^ each other, leading to the same growth rate. 
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Figure 7. (a-c) Illustration of By (z, t) from non-rotating turbulence, using 
CTh = (Tu with S = 2, Lz = 16, for DNS, DQLS, and CE2 from (a)-(c). 
These parameters are identical to I^. ^ but with cri, = cTu = cr/2, where 
cr is the driving noise used in Fig.j^ Krt (d) shows the growth in time of 
the mean-field for (solid,blue) nonlinear equations; (dashed,red) quasi-linear 
DNS (the dotted black line shows the fluid energy in the nonlinear run), (e) 
mean-field growth from CE2. As shown by the CE2 growth, the coherent 
dynamo is much stronger than in the kinematic case, but because the direct 
simulations (a-b) start from high amplitudes, it is hard to see the exponential 
dynamo growth phase in these simulations. 


homogeneity in the fluctuations u and b, such that V x 6 9 ^ 0. 
Since we assume initially homogenous fluctuations (termed 
Mo and ^ 0 ), this inhomogeneity must be introduced by B, 
which is assumed small. Considering the linearized fluctu¬ 
ation equations for simplicity (this is just the quasi-linear dy¬ 
namo, which we know works in any case){^it is evident that 
the kinematic dynamo arises from inhomogeneity induced in 
b fluctuations through the term V x (mq x B) in the fluctuation 
induction equation [Eq. ( |7b| )] . This leads to an inhomogenous 
contribution to 6 through mq x Z^mhom- In contrast, in the pres¬ 
ence of b{), an inhomogenous part of u will arise from the 
Lorentz force b^ 'VB + B 'Vb^ [see Eq. (|7a|)], giving a con¬ 
tribution to 6 through Minhom X b{). Without a mean-fleld flow, 
such a contribution is not possible from the induction equa¬ 
tion a^ne. In calculating the transport coefficients (Secs. 4.1 
and |5.2| ) we have verified that artificial removal of the Lorentz 
force causes the transport coefficients to return to their kine¬ 
matic values. It may be interesting in future work to examine 


Inclusion of nonlinear terms intro duces several additional complexities, 
see |Rheinhardt & Brandenburg |2010| . 


Figure 8. Same as Fig. [S] but using Keplerian rotation (cf. Fig. [^. The 
behavior is broadly similar to Fig.[8] with a slightly higher coherent growth 
rate because the velocity fluctuations have a positive effect in this case. 


in the vorticity dynamo (i.e., generation of U) in more detail, 
in particular its inte raction with the magnetic dynamo ( |Cour- 


voisier et~n^|2010| ). These effects are almost certainly much 
more important m the non-rotating case. 

Before proceeding it is worth commenting on an impor¬ 
tant difference between the magnetic shear-current effect dis¬ 
cussed below and the standard magnetic a-effect. This differ¬ 
ence stems from the fact that the magnetic a-effect can have 
either sign, since it is related to the small-scale current helic- 
ity, aM ~ -{b 'V xb). In practice, as the small-scale dy¬ 
namo grows in the presence of helical velocity fluctuations. 
au grows with the opposite sign to the kinematic a-effecf^ 
the origin of catastrophic quenching ([BTackman & Eield 


2002} [Brandenburg & Subrarrianian|j2005 

the rnagnetic shear-current effect drives t 


. In contrast, since 

le dynamo through 


a resistivity, 17 ~ {b^), its sign is fixed. This implies that the 
the source of magnetic fluctuations can be the small-scale dy¬ 
namo, in some sense the inverse of qu enching. In a recent 
paper ( [Squire & Bhattacharjee||2015b| ), we have shown that 
this mechanism is realizable at higher Rm where the small- 
scale dynamo is unstable. In particular, we see a decrease in 


One possible exception to this that may be very important could occur 
in the presence of magnetic instabilities, for instance the MRI. In this case it 
seems more likely that the magnetic a effect might overwhelm the kinematic 
effect, since b fl uctuations do not arise purely as a conse quence of small-scale 
dynamo action ^Gressel|2010||Park & Blackman|2012^. 
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riyx after saturation of the small-scale dynamo, which can in 
turn drive a coherent large-scale dynamo. 

5.1. Numerical experiments on the magnetic dynamo 

Here, we argue for the existence of the magnetic dynamo 
through numerical experiments, in a similar way to the dis¬ 
cussion in Sec. Since there is no small-scale dynamo due 
to the low Reynolds numbers, we excite homogenous mag¬ 
netic fluctuations {b^) by forcing the induction equation with 
(Th (the statistical properties of which are chosen to be the 
same as 

Figures [7] and illustrate the comparison between DNS, 
DQLS, anaCE2, at identical parameters to Figs. and and 
using the same total forcing level (i.e., cTb = (Tu = (rjl, ^ere 
(T is the driving noise used in Fig s. [T] and [2|). The most obvi¬ 
ous difference - comparing Figs. Hand [^o Figs. |T] and |^- 
is the much higher amplitudes in the direct numerical simu¬ 
lations (both quasi-linear and nonlinear). This is not due to 
the mean-field dynamo and simply results from the approxi¬ 
mate equipartition of B with b due to the finite size of the do¬ 
main. This occurs almost immediately because of the strong 
b{). Thus, the strong magnetic fields observed at later times in 
the direct simulations are in the nonlinear saturation regime 
of the dynamo, where pyx might be expected to change sign. 
Since this paper is concerned with the linear growth phase, we 
shall not analyze this saturation phase in detail. 

In contrast, CE2, by eliminating all fluctuations in 6, al¬ 
lows the mean-field exponential growth phase to be observed 
despite the presence of strong magnetic fluctuations. (We re¬ 
mind the reader here that the fundamental model used in CE2 
is identical to DQLS, the only difference arises from statis¬ 
tics being directly inserted into 6 to drive the mean-fields.) 
Comparing with Figs.[TJd) andj^d), we see that in both cases 
the mean-field growth is substantially faster; that is, the mag¬ 
netic fluctuations are contributing significantly to mean-field 
growth through the shear-current effect. In fact, since the 
growth rate is still strong in the non-rotating case [Fig.jTJe)], 
it is clear that the magnetic pyx significantly overwhelms the 
(positive) kinematic pyx, for the same forcing level cTu = CTb- 
Thus we have a mean-field dynamo driven by the magnetic 
shear-current effect. Moreover, the magnetic fluctuations pro¬ 
duce a stronger dynamo driving than the velocity fluctuations. 

Quantitative comparison of the CE2 calculations with the 
direct simulations is not possible due to the importance of 
nonlinear effects on the dynamo in the direct runs. While there 
appears to be an exponential growth phase in each direct case 
(before t ^ 100), it is somewhat too short to say for sure. In 
both the nonlinear and quasi-linear runs the wavelength ap¬ 
pears to be significantly shorter (k^ 47r/L^ ^ 67r/L^) than 

in CE2, for which the growing mean-field is the largest mode 
in the box. Since a stochastic-a effect is expected to be im¬ 
portant (presumably at a similar level to the kinematic case), 
this is not surprising; a fluctuations will act to increase the 
growth rate, decreasing the wavelength of the most unstable 
mean-field mode. The nonlinear and quasi-linear evolutions 
are broadly similar, although the nonlinear runs saturate at 
slightly lower amplitudes than the quasi-linear cases (a de¬ 
tailed comparison is not possible without running an ensemble 
of such simulations). Given this similarity - combined with 
the knowledge that the DQLS mean-field is, at least partially, 
driven by a coherent effect (the CE2 dynamo is unstable) - we 
conclude that this magnetic shear-current effect should also be 
playing a significant role in both the rotating and nonrotating 
direct numerical simulations [Figs.jTfa) and[^a)]. 


NONLINEAR DNS Knowing that the magnetic fluctuations can 
drive the coherent mean-field dynamo, it is helpful to examine 
the qualitative changes that occur as bo is increased. As a 
simple numerical experiment with the nonlinear equations, we 
start from pure velocity forcing and increase the driving in the 
induction equation, while keeping the total forcing, (Tu+o-j, = 
o", fixed. (While we have carried out these experiments both 
with and without Keplerian rotation, we present only the non¬ 
rotating cases here as the rotating results are similar.) 

Results in the range (Tb = 0 ^ 0.5o- are illustrated in Fig.[^ 
At (Tb = 0, we see a similar dynamo to that in Fig. [J 
though it is a little weaker due to the lower Urms and cnoice 
S = This is a stochastic-a effect, as seen by the slowly 
growing mean-fields that wander significantly in phase. Let us 
now consider the more interesting behavior of the other cases, 
(Tb = O.lo", (Tb = 0.2(t and (Tb = 0.5(r. Firstly, note that the 
larger mean-fields compared to the kinematic case are purely 
due to equipartition of B with b, as in Fig. [^a). Instead, 
our main result is the substantial qualitative difference in the 
appearance of the mean-field evolution between (Tb = 0.1(r 
and the cases with higher magnetic forcing. Specifically, at 
(Tb = 0.1 (T one observes a wandering mean-field as well as 
possibly a slow growth, behavior we interpret as a stochastic- 
a effect near its saturated state. In contrast, at (Tb = 0.2(t and 
(Tb = 0.5o", a relatively fast growth of B is observed until sat¬ 
uration at substantially larger values than seen at (Tb = 0.1 (t. 
In addition, the profile of By (z, t) for (Tb > 0.2o" is relatively 
constant in phase, suggesting that the dynamo is coherent. 

This behavior again suggests that a coherent dynamo can 
be driven by small-scale magnetic fluctuations - the magnetic 
shear-current effect. This dynamo saturates at larger held 
strengths than the stochastic-of dynamo, with the saturation 
amplitude being roughly independent of the level of magnetic 
fluctuations [as seen by comparisorp^ of Figs.|^c) and (d)]. 
Note also that this dynamo held appears to show quasi-cyclic 
behavior of some sort in its nonlinear regime [in Fig.j^d) the 
large scale held reappears again at later times]. The reason 
for this interesting behavior and i ts relevance to other dyna mo 
cycles (e.g., in MRI turbulence, |Lesur & Qgilvie||2008b| ) re¬ 
mains unclear, and given its apparent origin m nonlinear dy¬ 
namo physics, we leave its study to future work. 

To ensure the observed behavior is robust, we have rerun 
each of the simulations in Fig. [^several times, varying the ini¬ 
tial conditions and random number seed. These (not shown) 
have illustrated that the (Tb = O.lo" simulations occasionally 
excite the coherent dynamo similar to that in Fig.[^c-d), and 
will eventually do so if evolved for a sufficiently long time. In 
addition, the (Tb = 0.2 (t occasionally fails to excite the coher¬ 
ent dynamo as quickly as observed in Fig. |^c). This brings 
us to the conclusion that the coherent dynamo can be excited 
for (Tb > (0.1 ^ 0.2) X (T and the simulation outcome de¬ 
pends on properties of an individual realization around this 
boundary. We have failed to find coherent dynamo excita¬ 
tion at (Tb = 0.05(t, having tested a number of realizations 
over very long time periods. This dependence on realization 
is very similar to the behavior observed in shear-dynamos at 
higher Rm, where the small-scale dynamo acts as the source 
of bo fluctuations ( Squire & Bhattacharjee|2015b| ). 


5.2. Direct calculation of transport coefficients 


Note that the different z domain size must be taken into account to com¬ 
pare Fig.with Fig.|^a). 
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Figure 9. Low Rm driven DNS at = 1, = 8, and no rotation, with cTu + crt = cr chosen to be constant in each simulation (the level is such that Urms ~ 0-2 

when only velocity forcing is used). From left to right we take (a) = 0, (b) cr^ = O.lcr, (c) = 0.2cr (d) = 0.5cr. The top row of each subfigure illustrates 

the time development of By (z, t), the bottom row illustrates the kinetic energy (dashed, blue) and magnetic energy (solid, red). 


As in Sec. |4.1| we can directly calculate the transport coeffi¬ 
cients of the magnetic dynamo by fixing the mean-fields and 
driving magnetic fiuctuations. Within quasi-linear theory, this 
is a straightforward generalization of kinematic calculations, 
and the transport coefficients in the presence of both magnetic 
and velocity fiuctuations will be the sum of those calculated 
with one or the other, rj = rju + rj^. However, inclusion of 
magnetic fiuctuati ons in the nonlinear test-field me thod can 
be more complex ( [Rheinhardt & Brandenburg|2010| ), and the 
linearity of the transport coefficients is lost, rf riu + lb (al¬ 
though of course at low Rm nonlinear results must approach 
the quasi-linear results). Because of this, we present only 
quasi-linear results for the magnetic dynamo coefficients, and 
leave magnetic test-field method studies to future work. 

Fig. [TOl illustrates rjyx and rixx when only magnetic fiuctua¬ 
tions are present, as calculated by setting au = ^ and fi xing 


By in CE2, with the same technique as detailed in Sec. 4.1 


Most notably, we see both that riyx and rixx are negative, both 
with and without rotation, and are of similar magnitudes. Im¬ 
portantly, a comparison of Fig. [T^ with Figs. [5^ shows that 
riyx is substantially larger in magnitude than the kinematic 
value, which implies that when brms ~ Urms the magnetic 
contribution should dominate. For example, without rota¬ 
tion, the quasi-linear magnetically driven rfyx is larger than 
the quasi-linear kinematic rjyx by approximately a factor of 2 
at Rm = 100, meaning the presence of magnetic fiuctuations 
could change the sign of rjyx and excite a coherent large-scale 
dynamo once brms ^ This prediction is not far off 

the observed transition at o-j, ~ 0.2(r in Fig. with the dis¬ 
crepancy presumably arising due to inaccuracies in the quasi- 
linear approximation, as well as the additional presence of an 
incoherent dynamo mechanism. Note that r]xx in the magnetic 
case is much smaller than the kinematic value and will cause 




Figure 10. Transport coefficients for the magnetic dynamo -rj^x (solid, blue) 
and -rjyx (dashed, red) as a function of Re (at Pm = 1), for ^ = 2 and 
(a) D = 0 (b) Keplerian rotation Q = 4/3. (Note that both r]xx and rjyx 
are negative). The calculations are carried out at = 4 using CE2 (as for 
the kinematic case, there is very little draendence of L^). Coefficients are 
normalized by the Urms values from Figs.j^andl^ such that the values of the 
u and b transport coefficients can be directly compared. 


only a very small (probably unnoticeable) change to rjxx un¬ 
less Urms ^ i^rms- This is basically in agreement with the 
well-known result that magnetic fiuctuations do not signif- 
















































12 


icantly quench the turbulent resistivity. [In analy tic SQCA 
calculations with the shear added perturbatively ( [Squire & 
|Bhattacharjee|201 5c| ) the contribution of brms to is exactly 
zero]. 


Overall, we see that results of Fig. 10 agree well with our 
conclusions from earlier in the section and from Figs. 00 
and[^ Magnetic fluctuations in the presence of shear cause a 
signmcant negative contribution to riyx, which can overwhelm 
(or enhance in the presence of rotation) the kinematic coeffi¬ 
cient. Thus, with sufficiently strong magnetic fluctuations, a 
non-helical coherent large-scale dynamo is possible through 
the magnetic shear-current effect. 

6. DISCUSSION AND CONCLUSIONS 

The main purpose of this work has been to propose and ex¬ 
plore numerically a novel possibility for large-scale magnetic 
field generation in turbulent plasmas - the magnetic shear- 
current effect. The basic idea is that in the presence of large 
scale velocity shear, small-scale magnetic fluctuations pro¬ 
duce an off-diagonal turbulent resistivity (jiyx) with the cor¬ 
rect sign to cause mean-field dynamo instability when coupled 
with the shear. This is the magnetic analogue of the contro¬ 
versia l shear-current effect ( [Rogachevskii & Kleeorin|[2003 


2004| ) and the O x / (or Radler) effect (in the presence of 


shear). Importantly, this effect opens the possibility of the sat 
urated small-scale dynamo driving the large-scale dynamo, in 
stark contrast to standard a-quenching ideas where the small- 
scale dynamo is harmful to mean-field growth. Reassuringly - 
and unlike the kinematic shear-current effect - the sign of the 
magnetic effect agrees between analytic SOCA calculat ions 
([Squire & Bhattacharjee||2015c| ), the r-approximatio n ^Ro- 
gachevskii & Kleeorin|2004[ ), and quasi-linear theory ( [Singh 
& Sridhar[[2011t magnetic results presented here). In addi¬ 


tion, all three closure methods agree that the magnetic effect 
is substantially larger than the kinematic effect (for similar 
fluctuation levels b rms ~ Urms)^ perturbat ive MRI shearing 
wave calculations ( [Lesur & Ogilvie[[2008a[ ) have also found 
similar results. We hope that this agreement speaks to the ro¬ 
bustness of the effect in comparison to its kinematic cousin, 
at both high and low Reynolds numbers. 

In addition to the magnetic dynamo, we have presented re- 


viously by 

a number of authors (e.g., [Rogachevskii & Klee- 

orinJ2003[ 

Yousef et al. 2008 al [Brandenburg et al. 2008 Singh 

& Jmgade 

[2013). Our primary result is the qualitative (and 


_ by the w ( 

n X / (or RMler) effect ( [Krause & Radler|1980[ ), which for 
anticyclonic rotation will cause the off-diagonal res istivity rfy^ 
to have the req uired sign for a mean-field dynamo ( [Moffatt &\ 
[Proctor[[l982[ ). We have seen in a variety of examples how 
this can cause a change in the mean-field dynamo from being 
completely driven by fluctuations in a (the stochastic-a ef¬ 
fect), to being at least partially driven by the off-diagonal tur¬ 
bulent resistivity. The change is observable both qualitatively, 
in the spatiotemporal evolution of By, and quantitatively, in an 
increase in the dynamo growth rate. 

This paper has focused on the d ynamo at low Reynold s 
numbers, similar to that studied by Yousef et al.[ ( [2008a|b[ ). 


This choice has the advantage of both removing the complied 
tions of small-scale dynamo from the prob lem, and enabling 
the use of the quasi-linear approximation ([Sridhar & Subra- 

' 5a[) with some de¬ 

advantage allows clean and 


i|2009t [Squire & Bhattachaij^[2015a[) 

t accuracy The former advantage all 


straightforward separation of kinematic and magnetic effects, 
while the latter enables the use of statistical simulation tech¬ 
niques (CE2) that make the differences between incoherent 
and coherent dynamos particularly transparent. Nonetheless, 
precisely by enabling these simplifications, the low Reynolds 
number case is also less interesting. In particular, the mag¬ 
netic fluctuations cannot arise self-consistently through the 
small-scale dynamo, which is far more natural than a direct 
forcing of the induction equation (except perhaps in the pres¬ 
ence of magnetic ins tabilities such as the MRI). To rectify 
this, in a recent paper ( [Squire & Bhattacharjee[2015'b| ) we give 
numerical results that illustrate that the magnkic fluctuations 
arising from the small-scale dynamo can indeed cause a co¬ 
herent large scale dynamo through riy^. 

Given the historical controversy surrounding some aspects 
of the shear dynamo, we feel it helpful to give a brief sur¬ 
vey of the relationship to several previous works. As men¬ 
tioned in the main text, our results here on the ki nematic 
dynamo agree v ery nicely with numerical results in Yousef[ 
et al.[ p008a|b[ ). In particular, our conclusion that rotation 


fu ndamentally changes th e shear dynamo is nicely supported 
by [Yousef et al.[ ( [2008a[ ) Fig. 5, and can even be observed 
in the spatiotemporal plots of their Fig. 4. We also find ba¬ 
sic agreemen t with the quantitative results of [Brandenburg[ 
et alj ( [2008[ ), for instance, the transport coefficient calcula- 


tions showing riyx > 0, since these are carried out kinemati 
cally (neglecting the Lorentz force). However, we tentatively 
propose a different interpretation of their Fig. 8 (and possi¬ 
bly Fig. 7), whereby the magnetic shear-current effect is act¬ 
ing to drive the observed mean-field dynamo coherently (note 
the high Pm, which should lead to strong magnetic fluctua¬ 
tions). In support of this we note the very coherent appear¬ 
ance of the dynamo, as well as the near cyclic behavior in 
the saturation phase (cf. Fig. [^. Of course, more work is 
needed to assess this possi bility more thoroughly. Similarly, 
the Rm > 1 simulations of [Singh & Jingade] ( 2013^ (Figs. 6- 
8) may permit a similar explanation, although it is unclear 
whether there is truly a small-scale dynamo here. Finally, 
we mention again the analytic work of [Rogachevskii & Klee 


orm 


( 2004[ ), where the magnetic shear-current effect is derived 


within the r-approximation, although the authors do not com¬ 
ment on the result extensively. Specifically, it is clear from 
their Fig. 3 that the magnetic effect is far stronger (when the 
mean-field is zero) than the kinemati c effect, in broad agree- 
ment with our results in this work and[Squire & Bhattacharjee[ 

dMBcl ). 

Of course, since this work has explored only the low 
Reynolds number regimes, a variety of future studies will 
be important. While we have illustrated that the magnetic 
shear-current effect can arise fro m the small-scale dynamo in 
Squire & Bhattacharjee[ ( [2015b[ ), more work will be needed 


to more precisely assess regimes in which the effect may pre 
vail. Of particular interest will be the interaction of the ef¬ 
fect with magnetic helicity conservation arguments. This has 
been explored analytical ly and using quenching mo dels in[Ro- 
;achevskii et al.[ ( [20Q6[ ) (see also the appendix of IBranden- 
>urg et al.[2008[), but more numerical studies would be needed 


manian[ __^ 

gree of accuracy 


before any definite conclusions can be drawn. It would also 
be interesting to explore the relevance of the magnetic shear- 
current effect in flows with helicity and a deterministic a ef¬ 
fect. Is it possible that the effect could be present, perhaps 
after saturation of the aQ. dynamo? This may also be com¬ 
plicated by recent results showing that shear may help to en¬ 
hance helical dynamos by reducing the small-scale field gen- 
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eration ( [Cattaneo & Tobias|2Q14|[T^ias & Cattaneo|2014| ). 

Finally, we note the likely applicability of the magnetic 
shear-current effect to self-sustaining magnetorotational tur¬ 
bulence, where magnetic fluctuations are often substantially 
stronger than velocity fluctuations. With the confluence of 
magnetic fluctuations and anti-cyclonic rotation, it seems 
reasonable to surmise that the magnetic shear-current effect 
should be important. Dynamo cycles observed in unstrati¬ 
fied magnetorotational turbulence bear some resemblance to 
the quasi-periodic behavior the saturated state of the magnetic 
dynamo (Fig. 0, and it has been concluded previously that 
the dyna mo arises through a negative rfy^ (|Lesur & Ogilvie 
|2008a|b| ). In addition, some of the most solid evidence tor the 
effect’s importance come s from the CE2 simulations in |Squire[ 
& Bhattacharjee|p015a|). Here, since the kinematic effect is 


far too weak and incoherent effects are excluded, the magnetic 
shear-current effect is the only possible mechanism to drive 
the dynamo. The agreement between the saturation of the dy¬ 
namo in CE2 and nonlinear self-sustaining MRI turbulence 
simulations, in particular through the Pm dependence, pro¬ 
vides solid evidence that the MRI dynamo is indeed driven 
by the magnetic shear-current dynamo studied in this work. 
Interactions betwe en this effect a nd the a-effect due to ver¬ 
tical stratification ( [Gressel||201Q[ ) may help to provide sim¬ 
ple mean-field models that could be helpful in observationally 
useful disk models. 

Whatever the outcome of the variety of questions proposed 
in the previous paragraphs, given the generic presence of ve¬ 
locity shear flows and magnetic fluctuations in astrophysical 
plasmas, it seems likely that the proposed effects should find 
application across a wide variety of objects and phenomena. 
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APPENDIX 


A. STOCHASTIC-Of SHEAR DYNAMOS: SOME NOTES ON PREVIOUSLY PROPOSED MECHANISMS 


There has been a wide variety of literature on stochastic-a dynamos in shear flows. Here we consider the relationship between a 
number of these works, and explain some fundamental differences that would have important consequences for their observation 
in simulations. We feel that this discussion is suitable for presentation in this work, since our primary purpose has been to propose 
an alternative to the stochastic-a mechanism. 

At least two fundamentally different dynamo mechanisms ar e possible from fluctuations in th e a effect with z e ro mean. The 


first, which has been explored for a variety of perspectives in [Vishniac & Brandenburg 


et al. 


McWiI rams 


instability. 


2008|);[^shby & Proctor] ( [MT0| ); ^einemann et^ 
(|2012|), has the property (discussed in Sec. 
II 


|20li|);|Richardson & Proctor 


We 


H]) that 


(fTWTl); 

2^; I 


Proctor] (2007]); 
jMitra & Brande 


Brandenburg] 

:nburg' p()121 ); 


wil 


that {B{t)) decays in time, and only {B^) undergoes exponential 
(We remind the reader that <•) refers to an ensemble 


term this the incoherent stochastic-a mechanism. 
average, while ~ refers to the me an-field average.) The second mechanism, which i s in essence the Kraichnan-Moffat dynamo 
(]Kraichnan]] 19761 ]Moffatt]]1978]), has been explored in the context of shear flows in|Silant’ev|(]20001);]Sridhar & Singh] (]2014]), 
and does exhibit growth of {B(t)). We shall term this the coherent stochastic-a mechanism. Since part of our argument for the 
prevalence of coherent dynamo in some of our numerical experiments has centered on the requirements on mean-field evolution 
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imposed by = 0, it seems worth explaining in more detail the coherent stochastic-of mechanism and its relation to the 

incoherent variety. __ 

In its absolute simplest form, the dynamo in |Kraichnan| ( [1976] ); |Sridhar & Singh| ( |2014| ) can be described as resulting from 


dtB = Vx [a(x, t)B] + r]TV^B, 


(Al) 


where a(x, t) is a spatiotemporal fluctuating a-effect, assumed to arise from smaller scale fluctuations, and rjr is the turbulent 
resistivity. One then specifies that {a) = 0, {a{x, t)a(x', t')) = 2Jl(x,x')D(t, f), and forms the equation for (B) 


d,{B} = Vx(VMX{B})^rjKV\B}, 


(A2) 


where rjK = -^(0), and Vm = f(^ dT(a(x, T)Va(x, 0)). For sufficiently strong a fiuctuations, instability arise s for (B), because 

Tjx bec omes negative. Note that for such an instability the smallest scales of the mean-field grow the fastest. [Sridhar & Singh| 
( 20141 ) give a variety of interesting extensions to this model, including the effects of non-zero a correlation time Tq,, and shear 
(which changes the dynamo only if Tq, ^ 0 ). 

Why is it that this dynamo seems to be mean-field in the true sense - {B) grows exponentially - while this is not true for 
the incoherent stochastic-a dynamo? This question is important for understanding the shear dynamo, since a dynamo arising 
though this coherent stochastic-of mechanism will have very different properties. While it seems that all previous treatments of 
this dynamo have considered a spatiotemporal fluctuating a coefficient, this is not the fundamental difference. In particular, if we 

simply assert that a(x, t) = Qf(0 the dy namo can still exist with rjx = rij - ck^/ 2, although Vm = 0. The answer to this question is 
given in|Mitra & Brandenburg|p012|) Sec. 3.3, where they examine the effects of mutual correlations between a coefficients. In 
particular (now considering specifically a horizontal mean-field average such that we have only a 2-D system), they find that in 
the presence of mutual correlations between a coefficients. 


{aij{t)aki{t')) = 'D%5{t-t'), 


the ensemble averaged mean-field {B) = {{Bx), {By}) satisfies the equation 


d,{B) 


-S + k\D^yl - D7x) -kHriT + DZ ■ 




{B). 


(A3) 


(A4) 


Evidently, from Eq. in the coherent stochastic-a mechanism, axxO) = ayy(t), while ayx(t) = axy(t) = 0. This implies 
Dyy = = Hyyy = (q(^)/ 2, while all other vanish. Thus, one obtains exactly the same instability from Eq. ( |A4| ), since 

rjT - can be negative. 

We thus see that the coherent stochastic-a mechanism requires the rather specific situation of strong diagonal a fiuct uations, 
but very weak off-diagonal a fiuctuations (since = (a^_^)/2, and similarly for axy). While the exact result Eq. ( |A4| ) is only 
valid for a with no spatial dependence, it seems almost certain that similar conclusions will hold if spatial variation is also 
included. Is it realistic for the correlation between axx and ayy to greatly exceed the fiuctuations in ayx and axy (their difference 
must also overcome tjt)? Possibly, for instance if the fiuctuations in aij arose purely from fiuctuations in small-scale helicity, 
but this situation strikes us as unlikely. In any cas e, it seems that more w ork, both numerical and analytic (e.g., inclusion of ayx 
and axy in the much more thorough calculations of|Sridhar & Singh|2014|), would be needed to thoroughly assess this possibility. 
Overall, the confluence of factors against the coherent stochastic-a dynamo - the requirement for very strong a fiuctuations, the 
significantly adverse effect of off-diagonal a, and the fact that one would observe a mean-field that grows much faster on the 
smallest scales - leads us to conclude that this mechanism has probably not been observed in previous numerical experiments on 
shear dynamo. 






















